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. Abstract 

In this paper, motivated by applications in computer graphics and animation, we study the 
numerical methods for checking C'^— regularity of vector multivariate subdivision schemes with 
dilation 21. These numerical methods arise from the joint spectral radius and restricted spec- 
tral radius approaches, which were shown in [4] to characterize WjJ'— regularity of subdivision in 
terms of the same quantity. Namely, the (fe,p)— joint spectral radius and the (fc,p)— restricted 
, spectral radius are equal. We show that the corresponding numerical methods in the univariate 

■ scalar and vector cases even yield the same upper estimate for the (fc, oo)— joint spectral radius 

for a certain choice of a matrix norm. The difference between the two approaches becomes 
apparent in the multivariate case and we confirm that they indeed offer different numerical 

■ schemes for estimating the regularity of subdivision. We illustrate our results with several 

examples. 
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^ ■ 1 Introduction 

l> 

■ Subdivision schemes are recursive algorithms that starting with a coarser mesh in W determine 

\ the coordinates of the finer vertices c'''^^^ by local averages of the coarser ones 

C^; c('-+i) =5acM, r>0. (1) 

The subdivision operator Sa is linear and describes the local averaging rules. The locality of 
_ the subdivision and algorithmic simplicity of the subdivision recursion ensure that ([1]) is fast, 

, efficient, and easy to implement. These features make subdivision popular in computer graphics 

\ and animation, see [31 [71 [131 [5S] and the references therein. 

We restrict our study to the shift-invariant setting and study the multivariate vector schemes 
defined by the subdivision operator 

5acW - ^ A(.-2/3)c«(/3), r>0, 

mapping the space ^"(Z**) of vector-sequences indexed by Z'' into itself. The associated subdivision 
mask A = [A{a.))^^^^ j^-^s, N G N, is a finitely supported matrix- valued sequence. 

A challenging task is to provide a characterization of the regularity of the limits of subdivision 
recursion and to develop numerical methods for checking their regularity. Two prominent methods 
that characterize the W^-regularity, 1 < p < oo, k € N, oi subdivision in the shift- invariant setting 
are the so-called joint spectral radius (JSR) approach [SI [H] and restricted spectral radius (RSR) 
[HIS] approach. Generally speaking, the joint spectral radius approach characterizes the regularity 
of subdivision limits in terms of the joint spectral radius of a finite set 

A^{A\v, : ee{0,in 

of square matrices derived from the mask A and restricted to a common finite dimensional 
invariant subspace Vk- The sequences in Vk annihilate certain polynomial eigensequences of Sa 
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and, thus, the process of restricting of to Vk mimics the computation of discrete derivatives 
V''. The essence of the restricted spectral radius approach can be formulated as follows: if the 
operator Sa admits the factorization 



then the restricted spectral properties of the associated difference subdivision schemes Sb^. char- 
acterize the regularity of the underlying subdivision. It had been believed until recently that the 
JSR and RSR approachers are intrinsically different. In [HE], we showed that these two ap- 
proaches characterize the W^— regularity, 1 < p < oo, of subdivision in terms of the same quantity 
(fc,p)— JSR and differ by the numerical methods they yield for its estimation. Motivated by appli- 
cations in computer graphics and animation we compare in this paper only the numerical methods 
for approximation of (A;,oo)~JSR. The problem of computing the (fc,oo)— JSR is NP-hard, i.e, 
there is no polynomial-time algorithm for its approximation, see The brute force Branch- 

and-Bound algorithm dating back to [TU] is available. Its approximation of (fc, oo)— JSR is based 
on the estimate 



r € N, and allows us to approximate {k, oo)— JSR with an arbitrary precision. The estimate in ([3]) 
is independent of the choice of the matrix norm || • ||. The Branch-and-Bound algorithm is compu- 
tationally expensive and its rate of convergence is not known. See [TCI IT5] for a successful attempt 
to reduce the computational complexity of the Branch-and-Bound algorithm. The restricted spec- 
tral radius approach leads to linear programming, see [4 , that due to the usual properties of the 
spectral radii yields another upper bound for {k, oo)— JSR. The computations in [24 show that the 
lower bound in ^ is not always reliable. 

The main goal of this paper is to compare the numerical methods arising from the joint and 
restricted spectral radius approaches for approximation the (fc,oo)— JSR. In section [5T] we show 
that if we choose the matrix norm || • || in to be the infinity matrix norm, then the upper bound 
in ([2]) and the one obtained using the optimization approach in P| Section 4.6] coincide in the 
univariate scalar and vector cases. Note that the flexibility of the joint spectral radius approach 
is that the estimate in ([3]) is independent of the choice of the matrix norm || • |j. This in some 
cases leads to a sharper upper estimate for (fc,oo)— JSR, see example 13.51 or even allows for exact 
computations of the (fc, oo)— JSR, see [TBI [18]. The results of section [3T] also show how one can 
easily determine the entries of the matrices A^^- jy^ from the entries of the corresponding difference 
maks Bk- In section [3.21 we show that in the multivariate scalar or vector cases the Branch- 
and-Bound algorithm and the optimization method presented in [51 Section 4.6] are intrinsically 
different and yield different upper bounds for (fc, oo)— JSR regardless of the choice of the matrix 
norm in In Section [?751 we compare the method for estimation of the (fc, oo)— JSR from [55] 
and the one based on the properties of the difference masks on an example of a divergent scheme. 
To summarize, examples illustrate that it is impossible to prefer one of the numerical methods 
discussed in this paper over the other and they can be used according to one's personal preference. 

2 Notation and Background 

An element /i — (^i, . . . , /is) G is a multi-index whose length is given by \fi\ :— fii + ■ ■ ■ + jj,s 
and /i! := /zi! . . . /is!- For a — (ai, . . . , a^) G Z'* and /i = (/ii, . . . , /is) G Nq define 



Denote by e^, £ = 1, . . . , s, the ^— th standard unit vector of Mf and by Cj, j = 1, . . . , n, the j— th 
standard unit vector of R", respectively. Let ^"^'^ {"L") denote the linear space of all sequences of 
nxk real matrices indexed by Z''. In addition, let P^'^ (Z*) denote the Banach space of sequences 
of n X d real matrices indexed by 7/ with finite oo-norm defined as 



(2) 




(3) 



:= <^..a^^ 



l|C| 



OO 



sup |C*(a)| 



oo ' 



(4) 
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where |C(q;)|oo is the oo-operator norm on M" if c? > 1 and the oo- vector norm if c? = 1. Moreover, 
let ^^^''(Z") c be the space of finitely supported matrix valued sequences. Specific 

examples of such scalar and vector sequences is the scalar delta sequence S G ^ol^") and the vector 
sequence SeJ e £j''"(Z'') defined by 



S{a) 



1, a = 0, 

0, aeZ'\ {0} 



and SeJ (a) 



ej, a = 0, 
0, aeZ^\{0}. 



(5) 



Let A £ £q^"(Z*) be a finitely supported matrix sequence shifted so that suppA C [0,A^]* 
for some e N. The subdivision operator Sa ■ ^" (Z*) ^" (Z'*) associated with the mask 
A e ^^''"(Z") is defined by 



SAcia) = ^ A(a - 2/3)c(/3), a G Z" 



(6) 



The subdivision scheme then corresponds to a repeated application of Sa to an initial vector 
sequence c e i!" (Z'*) yielding 



The dimension of the following subspace of M" 



:= < w 6 K" : XI ^('^ ~ = w, e G {0, 1}" I 



(7) 



(8) 



determines the structure of the difference operators V'^' we define next. Let m := dim{£A)- For 
C G r ^''(Z^) and D G £'^^"(Z'*) we define the j-th column of V^C, : r^''(Z'*) r ^''(Z^) 
as 



1 < ^ < s, 1 <j <d, 



(9) 



and VfD :— (VeD^^ , respectively. Let fc G N. For our analysis we make use of the notion of 
the fc— th difference operator, the discrete analog of a derivative. The fc— th difference operator 
yk . (nxd^is-^ _^ fuN^.uxd^'^s^^ j^^^^^ ^ , is defined by 



V 2 . . . V g J„ 





In- 



(10) 



We say that the subdivision scheme 5a is — convergent, if for any starting sequence c G 
€JJ-,(Z'*) there exists a vector-valued function /c G {C^{W))'^ such that for any compactly supported 
stable test function g G C^{W) 



hm max p'7c - gin * {2^^^^^^' . . . S\c){2'- ■)\\^ = 0. 

r->oo |;i|<fc 



(11) 



For more details on the properties of test functions see [9]. 

Another useful tool for studying subdivision schemes is the Laurent polynomial formalism. For 
a finite matrix sequence A G ^q^'^ (Z*) we define the associated symbol as the Laurent polynomial 



A*(z):=2-^ X A(a)z", zG(C\{0}r 



(12) 



where, in the usual multi-index notation, 2" = 2"^ • 



3 



The concept of the joint spectral radius was first introduced in |28j and was rediscovered in 
|f 0| for characterizing the continuity of univariate subdivision hmits. For extensions of the concept 
of the JSR see [HI H^- Generally speaking, the joint spectral radius approach characterizes the 
Wp— regularity of subdivision limits in terms of joint spectral radius (JSR) of a finite set of 

linear operators. These are derived from the linear operators '■ ^^^"(Z*) — ^ ^g'*"(Z*) satisfying 

Aev = J2 ^ia)A{e + 2 • -a), v e ^o''"(^')' e ^ {0> 1}'' 

and then restricted to the invariant subspace Vk C ^^^"(Z'' ). The elements of the finite dimensional 
subspace 



Vk := {v e £i''"([0,iV]") : ^ t;(/3)w(-/3) = for all u e Uk}, 



(13) 



annihilate the polynomial sequences inUk G £"(Z*) reproduced by one step of subdivision recursion 
in ([7]). For details on the structure of Uk see j4]. Thus, the process of restricting to the subspace 
Vk imitates the computation of the partial derivatives of subdivision limits function, see [H Section 
2.3]. As the subspace Vk is finite dimensional we replace the finite set of operators 



A\v, -.^{AeW, ■■ ee{0,l}^} 

by a finite set of square matrices 

A, := [A^{e + 2a - /?)] ^^^^[o^^]. , e £ {0, l}^ 
and use them for estimations of the (k, oo)— joint spectral radius 

l/r 

r 

For a difference scheme Sski k > 1, satisfying ^ define the restricted {k,oo)-norm 



PooiAlvk) := lim max 



Bfc|v'=l|oo 



sup 



I ^B. Veil 
IIV^clloo 



ceC(Z^),V'=c^O 



and the restricted (fc, oo)— spectral radius ((fc, oo)-RSR) 

Po.{SbM:= lim \\SL\^4'J,\ 



(14) 
(15) 

(16) 

(17) 
(18) 



For the results on existence of Sb^ see e.g. |nill]. To estimate the (fc, oo)— RSR we use the standard 
properties of spectral radii 

P^iSsJ^) = inf IIS'bJvMI^^'' 

and get the following estimate 

Poo{SbM<\\S^bM]1'^ (19) 

for any r G N. The restricted norm HS'^^ lyfc ||oo can be computed using linear programming in [3J 
Section 4.6]. The main result of ^4, shows that (fc,p)— JSR and (fc,p)— RSR are equal. 



3 Methods for approximating (/c,oo)— JSR 

In this section we show how the entries of the matrices A;,, \ v^, and their products depend on the 

entries of the difference masks B), , r > 1 and 1 < fc < A^. In particular, we show that the upper 
bounds in ([3]) and (jl9l) on the (fc, cxd)— joint spectral radius coincide in scalar and vector univariate 
cases for the special choice of the matrix norm, i.e. 



'BfclV'lloo — ll-S'stlloo 



max 



n 



A, 



(20) 
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Varying the matrix norm in ([3]) can lead to a better or a worse upper estimates ol the joint spectral 
radius, see Example 13.51 In the scalar or vector multivariate case, the identity 



max 



does not hold in general, see subsection 13.21 Therefore, the JSR and RSR approaches lead to 
intrinsically different numerical methods for checking regularity of subdivision. 

An important observation that connects the entries of the matrices . | v;, with the ones of the 
difference masks Bk, k >1, states the following. 

Proposition 3.1. For k,r > 1, £j G {0, 1}'' and /3 £ Z'^ we have 

A. • . . V^J/„(- - ^) = ^ i?r(ei + • • ■ + 2'-i£, + 2'-/3 - ^)V'=5/„(- - 13). 

Proof. Let r > 1 and (3 E Z^. By [T71 Lemma 2.2] we have 

where a = ei + ■ ■ ■ + T'~^e\ + 2''7, 7 G Z''. Note that the application of Ae to Sl^bln means that 
the operator Ae is applied to each row sequence of V^SIn separately. Next using the identity ^ 
we obtain ^ ^ 

^B, V^-<5/„(a - /3) = . ■ . A. V'=<5/„(7 - /3), 7 e Z\ 

The definition of and the definition of the standard convolution operator * yield that the 

elements of the sequence . . . V'^'5/„(- — (3) have the following representation 



in terms of the elements of the subsequence B\, {si 
claim follows. 



2'-i£, 



2"- (7), 7ez^ 

+ 2'---/3) oiB^^\ Thus, the 

□ 



3.1 Univariate scalar and vector cases 



Let us see how to use the result of Proposition 13.11 to compare the upper bound in ([3]) and (I19p . 
For that we need the following auxiliary lemma. 

Lemma 3.2. For r > I, I < k < N , and S"^^ : ^"(Z) r (Z) we have 

II'S'bJv'=IIoo = II-SbJIoo- 

Proof. We consider the vector case first. Due to 12, Proposition 6.10], the non-restricted operator 
(fc, cxd)— norm \\S^^ \\oo is defined as usual by 



'Bt II 00 



and can be computed as follows 
\\S\ 



^Bk Hoc 



= max 

ae{0,2'--l} 



The absolute value of the matrices b''^^ in (|2T]) means that we take the absolute values of each 
matrix entries. Note that, due to the compact support of A and the definition of V'', we have 
s\vp^{B^^^) C [0,iV-|-A:] for all r > 1. Thus, ([21]) is equivalent to the following linear programming: 
for i = 1, . . . , n and a G {0, 2^ — 1} 



Y^\B^;\a-rp)\ 



1} 



(21) 



max Y.(^i%A<^-2^f3)dAf3) 

l3e[-N-k,0] j=i 



(22) 



-l<d,(/3)<l, /3eh7V-fc,0], j = l,. 
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For fixed i and a, the solution d*-*'"-* G £^{[~N — k, 0]) of this optimization problem is given by 

= sgn(i?i'^^)., {a - 2^/3) , Pe[-N-k,Q]. 

We get the maximizing sequence d G i^oii^^ ^ k,0]) for all i and a choosing d = d^*'"-* such that 
the number determine by linear programming in (j22p for this d is maximal. 

On the other hand, from 4^ the computation of the restricted (/c, oo)— norm is equivalent to 
the following linear programming: for i = 1, . . . ,n and a £ {0, 2*" — 1} 

n 

max ^ ^(i?('-)),,(«-2'-/3)V'=c,(/3) 

/3G[-Af-fc,0]i=l (23) 

-l<V^c,(/3)<l, /3g [-iV-/c,0], j^l,...,n. 

For fixed i and a, the solution c'^" e t^{[—N — 2k, 0]) of this optimization problem is given for 

c;^"(/3)=sgn(i?l:)),, 



each j by the linear system of equations V^c^^{l3) = sgn(i?^'^'')y (a — 2'"/3), j3 G [— iV — fc, 0]. The 



solutions Cj exist and are unique, due to the invertibility of the corresponding matrix of this system, 
which is bi-diagonal with —1 and 1 on the main and upper diagonals, respectively. We get the 
maximizing sequence c G P^{[-^N — 2k, 0]) for all i and a choosing c = c^*'") such that the number 
determine by linear programming in (|23p for this c is maximal. Therefore, p.2p is satisfied. 

In the scalar univariate case the proof is analogous, one should only take into account that 
the support of Bk decreases with k and replace above [0, N + k] by [0, N — k] and [—N — k, 0] by 
[—N + k,0], respectively. □ 

Combining Proposition 13. II and Lemma 13.21 we obtain the following result. 

Theorem 3.3. For r > 1 and 1 < k < N we have 

r 

e,?fc?i}lln^^>^IL=ii^kiv^L- (24) 

j=l 

Proof. We start with the scalar case. Due to ^ Lemma 4.4, 4.9], we have 

T4 =span{v'=<5(--/3) : /3 G [0, - fc] } . (25) 



Thus, due to the fact that the spanning set above is also a basis for 14 and by Proposition 13. 1[ we 
get for any ?- > 1 and ej G {0, 1} 



Therefore, 



Bl\ei + --- + 2'-^er + 2''p~ p) 



/3,/36[0,Ar-fc] 



£^£{0,1} 



J = l 



and the claim follows by Lemma [ 

In the vector case, define as in [U Lemma 4.9] the spanning set of Vk to be 

I (V"<5e,y (■-/?) : pe[Q,N-k], j = 1, . . . , m| y {,5eJ(. - /?) : /3G[0,iV], j = m + 1, . . . , n} . 

This set is also a basis for Vk ■ The corresponding matrix representations of Ae \ vi. are obtain as 
follows. First consider the set Vk spanned by 

I (V"5e,)^ (■-/?) : pe[0,N], j = l,...,m^[j[SeJi-^p) : p e [0, N], j = m + 1, . . . ,n] . 
By Proposition 13. II we get for any r > 1 and Sj G {0, 1} 



Bl:\ei + --- + 2'-hr + 2^f3~P) 



/3,/36[0,W] 
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As the spanning set of Vk is not a basis for Vk we need to remove the rows and columns of ^ely^^. 
that correspond to 



V'^Se.f {--13) : p e [N - k + I, N], j ^ 1, . . . ,m} 



Then the claim follows as in the scalar case. 



□ 



Remark 3.4. (i) The size of the matrices Ag|v^ in (jS]) is of importance for numerical computations 
of {k, oo) — JSR and is determined by the dimension ofVk- For a method that allows us to reduce the 
sizes of these matrices see 1271. In the scalar case, the restriction of to Vk in (j25p is equivalent 
to removing trivial cycles of A*(z). 

(ii) Note that the choice of the matrix norm on the left hand- side of ()24|) is crucial. For any 
other matrix norm the statement of the Theorem I j'.gl is not true in general, see examples below. 
This illustrates one of the advantages of the numerical method arising from the JSR approach and 
is exploited in JTW, where the authors approximate the so-called extremal matrix norm that allows 
for exact computations of the (k, oo) — JSR. 



Let us next illustrate the result of Theorem 13.31 on some examples. 
Example 3.5. Consider the A— point scheme with the mask given by 



A*{z) 



1 



-1 + 9z^ + 16z^ + 9z 



4 z^ 



1 



) = — (1 + z) (-1 + 2z + 6z^ + 2z^ - z^) 



32 ' '32 
The corresponding subdivision scheme is ll^l- Due to 

= 1 and A*(-1) = 0, 

by ^ Corollary 3.9] we can compute the difference masks Bi and B2 satisfying the equivalent 
formulation of ^ 

{z-l)''A*iz)^BUz)iz^-l)^ fc = l,2, 

and given by 



Bl{z) 



32 



(1 + z) B*{z), B*{z) = _ (-1 + 2z + + - z^) 



By Lemma 2.4] we have 



Vi = span|vJ(- - /3) : < /3 < sj 



and 



V2 = span|v2(5(- - /3) : < ^ < 4|. 

The corresponding matrix representations of Ae\vkJ ^—^t^, are 



Bi(0 + 2/3-;3) 



Ai\v, = + 2/3-/3) 



/3,/3e{0,5} 



1 

16 



1 

16 



-1 


8 


1 








" 





1 


8 


-1 











-1 


8 


1 

















1 


8 


-1 


1 


8 


-1 








' 


-1 


8 


1 














1 


8 


-1 














1 


8 


-1 






and 



Ao\v, = 52(0 + 2/3-/3) 



/3,/3e{0,4} 



1 

16 



0-1 







~1 

6 -1 



7 



^ll 



^2(1 + 2/? 



/J,/3e{0,4} 



1 

16 



2 2 
-16-100 

2 2 

2 2 



Hence, on the one hand we get 

\\SbA 



■[\\Aq\vA\oo. PilyJIoo} 



On the other hand 
and 



ll^'salloo = maxjllAolyJIoo, milvallooj = ^ 

ll^sJIoo <max{||Ao|yJ|2, IIAilyJb} =0.7195 
II^bJIoo > niax|||Ao|yJ|2, Pilvslb} =0 



= 0.467.. 



The second example is of a vector univariate subdivision scheme. 
Example 3.6. Let the mask A e (^'^'^{'L) he given by its symbol 



A*{z) = 



1/2 1/4 
1/4 



1 3/4 
1/2 



1/2 1/4 
1/4 



Due to 



A*(l)-(1,0)^ = (1,0)^ and A*(-l) • (1,0)^ = (0,0)^, 
by Q Corollary 3.9] there exist the difference mask Bi satisfying the equivalent formulation of ([2]) 



z-l 
1 



A*{z) = B*{z) 



and given by 



Bl{z) = 



1/2 -1/4 
1/4 



1/2 -1/2 1/2 
1/2 i ^ I 1/4 



Note that in the vector case, due to the definition of the difference operator V, the support of the 
difference mask Bi is larger than that of A. To compare the entries of Ag^lv^ to the ones of Bi, 
we first, as in the proof of Theorem \3.3\. define Vi to be the set spanned by 



{(V<5ei)^ (•-/?) : /3e [0,2]}u{5e^(--/3) : /3e[0,2]} 



Then we get 



My, 



Bi(£ + 2/3-/3) 



^i,/3G[0,2] 



ee {0,1} 



from which we obtain 5x5 square matrices A^jy^ , e £ {0,1}, by removing one but last rows and one 
but last columns of the corresponding matrices ^ely^. Note that such a choice of Ag\v^, s £ {0, 1}, 
depends on the structure of the basis of Vi give by 



{(V<5ei)^(--/3) : /3 G [0, 1] } U { Je^(- - /3) : /3e[0,2]} 



And we finally get 



\\Sbi\\oo = max|||Ao|yi lloo, ||^i|yi||oo| 
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3.2 Multivariate scalar and vector cases 



The result of Theorem 13.31 does not hold in general in the multivariate case. One of the reasons 
for that is that the spanning set of Vk in 4, Lemma 4.9] is not a basis for Vk, although the 
representation from Proposition 13.11 is still valid. Another reason is that the result of Lemma 13.21 
does not hold in general in the multivariate case, see |3J Example 5.2]. The following example 
illustrates that the numerical methods arising from the RSR and JSR approaches are indeed two 
different methods for estimating the (fc, cxo)— JSR. 

Example 3.7. Consider the scalar multivariate scheme Sa '■ ^(Z^) given by its symbol 

A*{zi, Z2) = ^ (l + ^i) (l + ^2) (1 + ^1^2) ■ 

Due to 

= 1 and A*{e)^l, e £ {1, -1}^ \ {1}, 
the first difference scheme satisfying the equivalent formulation of ^ 



zi-l 
Z2-I 



exists and is given, e.g. by the mask 



A*{z) = B*{z) 



1 




1 



1 
1 



1 

1 



1 




1 



As all the entries of Bi are non-negative, we apply ^ Remark 4.13] and get 



max IIS'bVcII 

l|Vc|U = l" 



max IIS'rcI 

|c|U = l" 



1/2. 



Note that, using \15\ Proposition 2.9], we can even compute Poo(S'biIv)- To do that we view Sbi 
as a scheme consisting of two scalar first difference schemes S^^ and S~ given by the masks 



Bi 





1/2 

1/2 1/2 

1/2 





with the entry at bold at the position (0, 0) 



Bi 







1/2 1/2 

1/2 1/2 





Then by '15, Proposition 2.9] we get 



Poo{Sbi Iv) = max |p {\Bi{a - 2^)]a,^gn) , P \ \Bi{a - 2fi)]o,^i3en 
where = [—2, 2]^ n 1? by ]15, p. 34-5] is any associated good set. A simple computation yields 



Poo{Sbi Iv) — 



1 



max IJS'bcI] 

c|loc=l 



Let us compare Poo(iS'b^]v) with the upper estimate for {0, 00) — JSR produced by ([3]). We fix a 
basis of Vi to be the following set 



{Vi5(--/3) : /3e{0,l}x{0,2}}U{V2^(--/3) : /3 = (1, 0), (1, 1)}. 
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Note that using Proposition \3.1\ we get 



Ai,o)V25(--(l,0)) = iv25. 



and 



(0,1)) + 1V2(5(.-(1,0)), 



2 2 2 

where V2i5 is not one of the above basis elements ofVi. Thus, the corresponding row o/^(i^o)|yi 
does not consist of the entries of the difference mask Bi as it happens in the univariate case. For 
all other possible choices of a basis of Vi we get analogous structure for at leats one of the rows of 
some of Ae\v-i- Therefore, we get 



max 



and also 



max 

e,G{0,iP 



3 

2> 



\Vi 



\S_ 



Bi Woo 



1.188. 



\SbAv\ 



JSR o 



> \\s 



Si IVlloo- 



Possible extensions of the result in |27j to the multivariate case is currently under investigation. 
Such an extension will not only allow us to reduce the size of A^^ \ y^. , but also may lead to different 
Vfc whose structure allows for better comparison of the RSR and JSR approaches. 



3.3 Divergence of subdivision schemes 

This scalar bivariate example is taken from \V7\. The mask is given by its symbol 

^ 1 /^l 3 2 3 1 2 

A {ZI,Z2) ^ - \-+Zi + -Z^ + -Z2 + Z1Z2 + -Z^Z2 

Note that the mask satisfies A*{1, 1) = 1 and A*{1, -1) = A*(-l, 1) = A*{-1, -1) = 0. It has 
been shown in [T| that 



PooiA\\ 



sup{|A|i/'-: r>0, X e a {A,,\v, ■ ■ ■ A,Jv^) , £3 ^ {OAV} 



(26) 



where a{M) denotes the spectrum, i.e., the set of all eigenvalues, of the matrix Af . In our case, 
one of the eigenvalues of y4(Q_o)|vi is equal to one, see [T7J page 1192] for details, and we get 



PooiA\ 



Vi 



> 1. 



This shows that the scheme is not uniformly convergent [5]. To get the same conclusion using the 
restricted spectral radius approach, we need to show that there exists i? S N such that for all r > i? 
we have ||5'bJv||oo > 1, which is a tedious task. Another possibility to show that the scheme is 
not uniformly convergent is to show that Sb^ does not converge to a zero limit function. To do 
that determine first the difference scheme Sb^ satisfying 



zi - 1 

^2-1 



A*{zi,Z2) = B*{zi,Z2) 



with the matrix-valued symbol B*{zi,zi) — 2 '^i^h the entries 



Kl{zi,Z2) 
^12(^11^2) 
b*22{zi,Z2) 



1 

16 
0, 
1 

16 



(1 + 3zi + 3z2 + Z1Z2), 

b*2^{z^,Z2)^^{2z2-2) 



16 



(z2 + 4zi+3). 



The symbol satisfies 



2 
2 
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and does not indicate that the associated scheme is not zero convergent. The problem is though 
that not all of the eigenvalues of the sub-symbols 

B*Jzi,Z2)^ ^Bi(e-2a)z", e G {0, 1}', 
at (1, 1) are less than 1, e.g 

Si,(o,o)(l,l) = ^ _2 4 ■ 

This violates a necessary condition for the convergence of subdivision, see [5]. Thus, the spectral 
properties of the difference scheme Sbi also indicates that the scheme is not uniformly convergent. 
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